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Beschreibung 

Anwendungsgebiet 

Die mechanischen Eigenschaften von biologischem Ge- 
webe (z. Elastizitatsparameter) sind fur die Beurteihing 
dcs Zustandes des Gewebes von groBem Ititeresse, In der 
medizinischen Diagnostik deuten Veranderungen der Hasli- 
zitatseigenschaften auf histologische und u. U, pathologi- 
sche Veranderungen hin. ALLgemein bekannt sind Prozesse 
wie die Bildung von Geschwiilsten und Verhartungen 
C'KnoteiT), die %, B. manuell tastbar sind ("Palpation"). In 
der Landwirtscbaft ist fur die Beurteilung der Qualitat des 
Fleisches von Schlachtvieb die Kenntnis mechanischer Ge- 
webeeigenschaften ebenfalls von Interesse. 

Stand der Technik 

Die sogenannte Tastbefundung ist ungenau und unemp- 
findlich. Wesentlich besser ist in dieser Beziehung die soge- 
nannte Elastographie, bei der elastische Gewebeeigenschaf- 
ten techoisch erfafit und z. B. in Form von Schnittbildern 
qualitativ oder quantitativ visualisiert werden, Dabei be- 
dient man sich hauplsachlich des Ultraschalls, wie er als 
bildgebendes Verfahren in der medizinischen Diagnostik an- 
gewandt wird, allerdings in einer rnodifizierten BetriebsarL 
In zeitlich nacheinander aufgenommenen Ultraschallbitdern 
konnen geringste Verschjebungen oder Verformungen in- 
nerhalb der dargestellten Gewebestruktur durch Auswertung 
der Bildsequenzen erfaBt und ausgewertet werden. Wird auf 
einen Gewebebereicb ein mechanischer Druck ausgeiibt, der 
eine Verformung des Gewebes 7.ur Folge hat, so verformen 
sich Bereiche mit unterschiedlichen Elastizitatseigenschaf- 
ten verschiedenartig. Das Elastographiesystem wertet diese 
Verformungen durch den numerischen Vergleich der Einzel- 
bilder aus und steiit die unterschiedliche Dehnungen im BiLd 
dar. Die notweadige Kompression des Gewebes wird extern 
durch den Ultraschallwandler provoziert oder intern durch 
Atmung oder Puis, Die Kompression ist nur gering und be- 
tragt bei Anwendung des ublichen diagnostischen Ultra- 
schalls nur Bruchteile von Millimetern; Wichtig dabei ist 
eine quantitative Kontrolle des Kompressionsdruckes, Der 
Druck wird in Richtung der Schallausbreitung provoziert 

Ein Verfahren der Ultraschall-Elastographie von Korper- 
gewebe ist erstmalig in einem Aufsatz von J. Ophir et al im 
Jahrc 1991 [1] T 12] bcschricbcn wordcn. Dabei werden Ul- 
traschaUbilder bzw. die korrespondierenden hochfrequenten 
Ultraschallechosignale, aus denen die ublichen Uitraschall- 
bilder in dem Gerat entstehen, so ausgewertet, daB "Verschie- 
bungen des Korpergewebes zwischen zwei, mit verschiede- 
ner Kompression aufgenomrnennen Gewebebildern berech- 
net werden, Auf diese Weise lassen sich die bereits erlauter- 
len Rucksehlusse auf die Elasiizitat des Organs bis hin zu ei- 
ner quantitativen Abbildung des Elastizitatsmoduls erzielen. 

Die hochfrequenten Echosignale aus einem komprimier- 
ten Gewebeabschnitt erreichen den Ultraschallwandler zu 
einem Zeitpunkt der von dem MaB der Kompression des be- 
treffenden komprimierten Gewebeabschnitts abhangt. Die 
Hauptaufgabe in der Auswertung der hochfrequenten Ultra- 
schallechosignale zur Berechnung der Dehnung besteht da- 
her darin, Zeitverschiebungen in kurzen zeitlichen Interval- 
len der hochfrequenten Ultraschalldaten zu berechnen. 
Diese Zeitverschiebungen zwischen korrespondierenden 
Echosignalen werden in alien durch das Ultraschallsystem 
abgebildeten Gewebebereichen in Form eines zweidimen- 
sionalen Bildes berechnet. Fur die Berechnung einer Zeit- 
verschiebung wird die Kreuzkorrelationsfunktion zwischen 
den hochfrequenten Echodaten verwendet, Einerseits kann 



die Zeitverschiebung durch Suclien des Maximums der 
Kreuzkorrelationsfunktion ermittelt werden [1]. Anderer- 
seits kann die Zeitverschiebung aus der Phase der Kreuzkor- 
relationsfunktion an irgendeiner Stelle, i. d, Regel an der 
5 Stelle Null, geschatzt werden [3] T [4] . 

Nach der Berechnung der Zeitverschiebung kann die lo- 
kale Dehnung durch Gradientenbildung oder mit Hilfe ein- 
facher li nearer Filter [6] berechnet und abgebildet werden. 



to 



Nachtelle des Standes der Technik 



Fur alle Verfahren zur Berechnung der Zeitverschiebun- 
gen wird die Kreuzkorrelationsfunktion korrespondierender 
Intervalle der Echosignale aus dem gleichen Gewebeab- 
15 schnitt unter verschiedener Kompression verwendet. Die 
Auswertung der Kreuzkorrelationsfunktion ist aufgrund der 
notwendigen Integration iiber das ganze Intervall zeitauf- 
wendig. Die bisher vorgeschlagenen Verfahren lassen sich 
in zwei Gruppen unterteilen: 

20 

1) Kreuzkorreiationsmethoden 

Verfahren, die das Maximum der Kreuzkorrelationsfunk- 
tion durch eine vollstandige Suche oder Einschacbtlung er- 

25 mitteln, Dazu miissen die Echosignale in der Regel mit einer 
sehr hohen Abtastrate abgetastet werden, um auch kleine 
Zeitverschiebungen genau messen zu konnen. Diese Verfah- 
ren sind auBerst zeitaufwendig, so daB sie nicht in ein Echt- 
zeitsystem oder Online-System implementiert. werden kon- 

30 nen, 

2) Phasenbasierte Verfahren 

Verfahren, die Zeitverschiebungen aus der Phase der 
35 Kreuzkorrelationsfunktion der komplexwertigen Uitra- 
schallechosignaie schatzen. Nachteiie dieser Verfahren sind 
mogliche Ungenauigkeiten und das Auftreten von Mehrdeu- 
tigkeiten bei groGen Zeitverschiebungen. Bisher konntcn 
diese Menrdeutigkeiten nur durch zusatzliche zeitintensi ve 
40 zweidimensionale Vor- oder Nachverarbeitungsschritte ge- 
wahrleistet werden. 

Aufgrund der Nachteiie der bisherigen Verfahren ist es 
nicht moglich, Ultraschaliechodaten hinreichend genau in 
Echzeit zu verarbeiten und somit Dehnungsbilder in Echzeii 
45 darzustellen. In der Vergangenheit hat sich aber gezeigt, dafi 
bei Ultraschallabbildungsvccfahrcn wic dem klassischcn B- 
Bild oder der Dopplersonographie die Echzeitfahigkeit der 
Systeme besonders wichtig fur eine breite Nutzung dieser 
Systeme ist 

50 Ein weiteres Problem ist, daB die Veranderungen der Si- 
gnallaufzeiten (Puls-Echo-Sequenz) die sich durch die 
Kompression des Gewebes ergeben, von der Entfernung 
zwischen dem Gewebeareal, weiche die Echos hervorrufen, 
und dem Ultraschallwandler abhangen, d. h. nicht konstant 

55 sind und mit. zunehmender Entfernung monoton ansteigen. 
Durch Schwankungen der Einhiillenden der Echosignale 
kommt es daher zu Ungenauigkeiten in der Verschiebungs- 
berechnung, da diese ein Interval] s bestirnmter Lange be- 
nutzL Diese wurden in der Vergangenheit u. a. durch Log- 

60 arithrnieren der hochfrequenten reellen Echosignale [7] ver- 
mindert Nachteilig wirkt. sich hierbei aber die Veranderung 
der Phasen der Signale aus. 



65 



Aufgabe der Erftndung 

Aufgabe der Erfindung ist es dahen eine Anweisung zum 
Handeln zu entwickeln, nut dem Ultraschaliechodaten hin- 
reichend genau moglichst schneil (z. B. in Echtzeit) ausge- 
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wertet werden konnen. Das Verfahren sollte die gleiche Ge- 
nauigkeit wie die Kreuzkorrelationsverfahren durch eine 
moglichst rechenefnziente Signalverarbeitung erreichen. 

Losung der Aufgabe 

In phasenbasierten Verfahren wird die Phasenverschie- 
bung zwischen zwei Echosignalen ais Ma8 fur die Zeitver- 
schiebung benutzt. Fur monofrequente zeitverschobene Si- 
gnale ist die Phase der Kreuzkorrelationsfunktion eine li- 
neare Funktion der zeitiichen Echoverschiebung. Bei der 
tatsachlichen, d, h. gesuchten Zeitverschiebung hat die 
Phase der Kreuzkorrelationsfunktion eine Nullstelle. Die 
Steigung entspricht der Kreisfrequenz der ntonofrequenten 
Signale. Dieser Zusamrnenhang kann fur die Berechnung 
der Zeitverschiebung aus der Phase der komplexen Kreuz- 
korrelationsfunktion an einer beliebigen Stelle genutzt wer- 
den. Die Phase der Kreuzkorrelationsfunktion kann aller- 
dings aufgrund der Mehrdeutigkeit der Arcustangensfunk- 
tion numerisch nicht eindeutig ermittelt werden, AuBerdem 
werden in Ulttaschall-Puls-Echosystenien moglichst breit- 
bandige Signale verwendet Fur bandbegrenzte Signale gilt 
obiger Zusamrnenhang nur approximative Diese Tatsachen 
fuhren zu den genannten Probiemen der Ungenauigkeit und 
der Mehrdeutigkeit phasenbasierter Verfahren. 

Beide Probleme (Mehrdeutigkeit, Ungenauigkeit bei 
Breitbandsignalen) konnen gelost werden, indem iterativ die 
Phase der Kreuzkorrelationsfunktion fur die Suche des Ma- 
ximums der Kreuzkorrelationsfunktion verwendet wird. Die 
vom Ultraschallwandler aufgenommenen Echodaien sind 
gemaB Abb, 1 in A-Linien organisiert. Die Berechnung von 
Zeitverschiebungen erfolgt an verschiedenen diskreten 
Punkten im Gewebe. In der Elastographie ist die Zeitver- 
schiebung, die an diesen diskreten Punkten berechnet wer- 
den konnen ein stetige, monotone Funktion uber der Hefe. 

Die Zeitverschiebungen werden sukzessiv fur jede A-Li- 
nie beginnend vom dem Ultraschallwandler nachstgelegen- 
den Punkt bis zum tiefsten Punkt berechnet, Fiir die Berech- 
nung der Zeitverschiebungen der Echos eines Punktes wird 
die folgende Iteration durchlaufen: 

1. ) Die Iteration wird begonnen bei der geschatzten 
Zeitverschiebung T, die fur den in der gleichen A-Linie 
dem Ultraschallwandler naherliegenden Punkt ermit- 
telt wurde. Fiir die Berechnung der dem UltraschabV 
wandlcr nachstgclcgcncn Punktc wird die Iteration bei 
Null begonnen, . 

2, ) Die Phase der Kreuzkorrelationsfunktion der kor- 
respondierenden Ultraschallechodaten an der geschatz- 
ten Zeitverschiebung T wird ennittelL Die ermittelte 
Phase ist aufgrund obiger Erlauterungen ungefahr pro- 
portional zurn Abstand der geschatzten Zeitverschie- 
bung zur wirklichen Zeitverschiebung. Die Proportio- 
nalitatskonstante entspricht in etwa der Schwerpunkts- 
frequenz des verwendeten Ultraschallwandlers (Kreis- 55 
frequenz). Die Schatzung der Zeitverschiebung X kann 
daher verbessert werden, indem der verbleibende Pha- 
senunterschied durch Teilen durch die Schwerpunkts- 
frequenz des Ultraschallwandlers in eine Zeitverschie- 
bung umgerechnet und die bisherige Schatzung der 60 
Zeitverschiebung entsprechend verandert wird. 

Durch Punkt 1.) wird gewahrleistet, daft die Kreuzkorre- 
lationsfunktion stets in der Nahe der wirklichen Zeitver- 
schiebung ausgewertet wird. So kann gewahrieistet werden, 65 
da8 verbleibende Phasenunterschiede eindeutig in Zeitver- 
schiebungen umgerechnet werden konnen. 

Durch eine mehrniatiges iieratives WiederhoLen von 2.) 



wird die Ungenauigkeit bei der Umrechnung von Phasenun- 
terschieden in Zeitverschiebungen bei der Verwendung 
breitbandiger Signale iterativ minimiert, Insgesami muB mit 
diesem Verfahren die Zeitverschiebung die Kreuzkorrelati- 
onsfunktion lediglich wenige Male (weniger als 10 Mai) fur 
die hinreichend genaue Schatzung einer Zeitverschiebung 
ausgewertet werden. Das vorgestellte Verfahren ist somit in 
zwei Beziehungen iterativ: Zeitverschiebungen werden 
durch wiederholtes Verwenden von berechneten Zeitver- 
schiebungen gleicher oder vorangehender Bildpunkte be- 
rechnet Da in der Elastographie auch sehr kleine Zeitver- 
schiebungen geschatzt werden miissen, die in der Regeln 
unterhalb der Abtastzeit liegen, miissen die Echosignale ge- 
eignet interpoliert werden. Bei dem vorliegenden Verfahren 
werden wie in [3] und [4] fur die Berechnung der Kreuzkor- 
relationsfunktion daher die komplexen, heruntergemischten 
Basisbandsignale verwendet. Da diese aufgrund der Band- 
paBbegrenzung von Ultraschall-Echosignalen eine deutlich 
geringere obere Grenz&equenz aufweisen, reicht hier eine 
rechenefnziente lineare Interpolation aus [5]« 

Ungenauigkeiten der geschatzten Phasenverschiebungen 
werden in diesem Verfahren dadurch vermindert, daJ3 die 
Amplituden der komplexen Basisbandsignale logarithmiert 
werden, nachdem sie zuvor urn eins erhoht wurden. Die 
Phasen werden nicht verandert. Dadurch kann die logarith- 
mische Kompression auch fur phasenbasierte Verfahren be- 
nutzt werden. 

Vorteile der Erfindung 

Das Verfahren erreicht im Gegensatz zu bisherigen pha- 
senbasierten Verfahren die gleiche Genauigkeit wie Kreuz- 
korrelationsmethoden. Dazu muB die Kreuzkorrelations- 
funktion nur einige Male (weniger als 10 Mai) ausgewertet 
werden, Auf diese Weise ist eine Echtzeitauswertung der 
Echosignale zur Berechnung der Gewebedehnung bereits 
jetzt mit geeigneten digitalen Signalverarbeitungsanlagen, 
wie sie in modernen digitalen Ultraschallgeraten zur VerfuV 
gung stehen oder separat als kommerziell erhaltlich sind, 
mdgHcfcu 
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Zeichnungen dargestellt und werden irn folgenden naher be- 
schreiben: 

Zu Abb, 2) 

5 

Schematische Darstellung eines System zur Aufnahme, 
schnellen Berechnung und von Darstellung von Dehnungs- 
bildern 

Das System zur Aumahme, schnellen Berechnung und 10 
Darstellung von Dehnungsbildern besteht aus drei Funkti- 
onseinheiten: 

a) ein handelsiibliches Uitraschallsystem mit verschie- 
denen Sonden. Diese konnen Ultraschallechodaten mit 15 
verschiedenen Scanprinzipien aufnehrnen (Sektor- 
Scan, Parallel-Scan). Ebenso spielt die Art der Fokus- 
sierung (eiektronisch oder durch einen Einzelwandler 
mit konstantem Fokus) keine Rolie. Die hochfrequen- 
ten Echodaten mussen in digitalisierter Form vorlie- 20 
gen. Bei Systemen mit analogem Beamforraer, oder bei 
Systemen mit Einzelwandler mussen die analogen 
hochfrequenten Echodaten direkt nacb dem Beamfor- 
mer, bzw. nach dem Einzelwander abgegriffen und di- 
gitalisiert werden. Bei Systemen mit digitalem Beam- 25 
former konnen die digitalen Echodaten direkt abgegrif- 
fen werden. 

b) eine Recheneinheit zur schnellen Berechnung der 
Dehnung. Dies ist das Kernstiick der Erfindung. Diese 
Recheneinheit kann sich innerhalb oder auGerbalb des 30 
Ultraschallgerates befinden. Die Recheneinheit kann 
dutch eine digitale SchaUung oder aber durch eine 
DSP-Karte, einen anderen Rechner., eine Recheneinheit 
des Ultraschallgerates realisiert werden (siehe Abb, 3) 

c) eine Darstellungseinheit. Diese Darstellungseinheit 35 
kann sich ebenfalls innerhalb oder auBerhalb des Ultra- 
schallgerates befinden. Es konnen die Darsteilungs- 
moglichkeiten des Ultraschallgerates verwendet wer- 
den. Bei AnschluB eines externen Rechner als Rechen- 
einheit konnen grafische Darstellungsmoglichkeiten 40 
dieses Rechners verwendet werden. 

Die Berechnung der Dehnungsbilder kann in Echtzeif, 
online oder offline nach der Aufnalime der Ultraschallcho- 
daten erfolgen. Fur die Darstellung der Dehnungsbilder 45 
{cann das Dchnungsbiid komplctt dargestellt werden, cin 
Teil des B-Bildes in einem besdmniteo Bereich (Region of 
Interest) durch das Dehnungsbilder iiberblendet oder einge- 
farbt werden. 

50 

zu Abb, 3 

ReaHsierung der RecbeneinheiL 

Die Recheneinheit kann in folgender Weise z. B r als Digi- 55 
talschaltung realisiert werden: zuriachst mussen die hochfre- 
quenten Echodaten quadraturdemoduliert werden (1). Die 
resultierenden Daten sind abwartsgemischt Sie werden Ba- 
sisbanddaten oder IQ-Daten genannt. Die Amplituden der 
IQ-Daten werden nun wahlweise logarithmisch kompri- 60 
miert (2), 

Die entsprechend aufbereiteten komplexen Echosignale 
(IQ-Daten) werden in einem Zwischenspeicher (3) abge- 
speichert, Dieser Zwischenspeicher muB die IQ-Daten je- 
weils zweier zeitlich aufeinanderfolgender B-Bilder spei- 65 
chern konnen. Er kann daher nach dem FTFOPrinzip im- 
mer, wenn die Echodaten eines neuen Bildes am Eingang 
anliegen, die Echodaten des vorletzten Bildes loschen. 
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Die Berechnung der Zeitverschiebung erfolgt in einer ge- 
sonderten Einheit (11). Diese kann in einer digitalen Schal- 
tung mehrfach vorhanden sein. Dadurch konnen die Zeitver- 
schiebungen mehrerer A-Linien gleichzeitig berechnet wer- 
den. Der Ablauf wird von einer globalen Ablaufsteuerung 
(12) gesteuerl. Diese Ablaufsteuerung steuert die Berech- 
nung der Zeitverschiebung an verschiedenen diskreten Posi- 
tioned Die Zeitverschiebungen werden gemafi Abb. 1 an 
verschiedenen diskreten Bildpunkten berechnet Verschie- 
dene A-Linien konnen parallel oder nacheinander abgear- 
beitet werden. Bei der Berechnung der Zeitverschiebungen 
einer A-Linie wird mit dem den Wandler am nachsten lie- 
genden Punkt begonnen und danach sukzessive Zeitver- 
schiebungen an defer liegenden Punkten berechnet Dieser 
Ablauf muB von der globalen Ablaufsteuerung (12) kontrol- 
liert werden. 

Die Einheit zur Berechnung einer Zeitverschiebung (11) 
enthalt folgende Elemente: 

Die komplexe Kreuzkorrelationsfunktion wird an der durch 
den Integrator (9) vorgegebenen Stelle ausgewertet (5). Der 
Integrator (9) dient als Speicher fiir die geschatzte Zeitver- 
schiebung wahrend der Iteration, Sein Startwert muB von 
der globalen Ablaufsteuerung gesetzt werden. Der Startwert. 
fur die dem UlTraschallwandler nachstgelegenen Punkte ist 
Null. Der Startwert fur alie anderen Punkte ist die Zeitver- 
schiebung, die fur den dem Ultraschallwandier nahergelege- 
nen Punkt berechnet wurde. Die Lage des Intervalls der Da- 
ten, die fur die Berechnung der Kreuzkorrelationsfunktion 
verwendet wird, wird von der globalen Ablaufsteuerung 
(12) bestimmt. Dieses ist ein Intervall konstanter Lange urn 
den Punkt an dem die Zeitverschiebung berechnet werden 
soli. 

Damit die Kreuzkorrelationsfunktion auch an Stellen be- 
rechnet werden kann, die nicht ganzzahlige Vielfache der 
Abtastzeit sind, mussen die Daten linear von einem Interpo- 
lator (5) interpoliert werden. Dieser greift auf die Basis- 
band-Echodaten der gleichen A-Linie der beiden B-Biider 
im Bildspeicher unter unterschiedlicher Kompression zu. In 
(5) wird somit die Kreuzkorrelationsfunktion der Basis- 
banddaten berechnet IJmdieKreuzkorreladonsfunktion der 
komplexen hochfrequente Daten zu erhalten, muB die Phase 
des Ergebnisses um das Produkt der Kreisfrequenz des Mi- 
schers (1) mit dem Zeitpunkt an dem die Kreuzkorrelations- 
funktion ausgewertet wird (9) in negative Richtung gedreht 
werden (6). Die Phase dieser komplexen Zafal wird nun be- 
rechnet (7) und mit dem rcziproken Wert der Schwcrpunkts- 
frequenz des Wandlers (Kreisfrequenz) multipliziert (8). 
Das Ergebnis wird zum Wert des Interpolators (9) hinzuad- 
diert Der Ablauf zwischen (5) und (9) rnuB einige Male 
(1-10) wiederholt werden. Dies wird von einer Ablaufsteue- 
rung (10) kontrolliert Das Ergebnis der Iteration wird an- 
schlieBend in einem Ergebnisbildspeicher (13) abgelegt Die 
Sleile, an der das Ergebnis dort abgelegt wird, wird von der 
globalen Ablaufsteuerung kontrolliert. Der Inhalt des Bild- 
speichers kann anschlieBend zur Anzeige gebracht werden. 

Samtliche Elemente dieser Recheneinheit konnen auf ver- 
schiedene Weise technisch realisiert werden: 

1. Jedes oder einige Element kann durch eine digitale 
Logikschaltung realisiert werden 

2. Jedes oder einige Elemente konnen von einer oder 
mehreren DSP-Schaltungen realisiert werden. 

3. Jedes oder einige Elemente konnen von einem 
Rechner innerhalb oder auBerhalb des Ultraschallgera- 
tes realisiert werden. 
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Patentanspriiche 

1 , Ein Verfahren und ein System zur Berechnung der 
Dehnung elastischen Gewebes mit. diagnostischem Ul- 
trascball, dadurch gekennzeichnet, daft die Zeitver- 5 
schiebung durch iterative Auswertung der Phasen der 
Kreuzkorrelationsfunktion zugehoriger Bchosignale 
geschatzt werden. 

2. Ein System nach Anspruch 1, dadurch gekennzeich- 
net, daB die iterative Auswertung der Phasen der to 
Kreuzkorrelationsfunktion an transduccrnahcn Gewe- 
bebereichen beginnt, bei denen die Nallstellensuche 
eindeutig ist und, daB hier gewonnene Ergebnisse bei 
der Berechnung von Verschiebungen bzw. Phasennull- 
durchgange in weiter entfernten Gewebebereichen in 15 
der Weise vorteilhaft genutzt werden konnen, daB die 
Kenntnis der vorher gewonnenen Phaseninformationen 
bei entfernten Bereichen vorhandene Uneindeutigkei- 
ten beseitigt. 

3, Ein System nach Anspruch 1, dadurch gekennzeich- 20 
net, daB die notwendigen Interpolationen hochfrequen- 
ter Bchosignale mit Hilfe einer linearen Interpolation 
von abwartsgemischten Echosignalen (Basisbandsi- 
gnalen) erfolgt. 

4. Ein System nach Anspruch 1 T dadurch gekennzeich- 25 
net, daB vor der die Genauigkeit der Berechnung von 
Zeitverschiebungen dadurch verbessert wird, daB die 
Amplitude der komplexen Bchosignale (unter Beibe- 
haltung der Phase) logarithmisch komprimiert werden > 
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